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The Anderson model with the Coulomb interaction between the local and conduction elec- 
trons is studied in the semiconducting phase. Based on a perturbation theory from the atomic 
limit, leading contributions for the c-f Coulomb interaction are incorporated as a vertex correc- 
tion to hybridization. An analytical solution shows that the effective attraction in the interme- 
diate states leads to a bound state localized at the local electron site. Self-consistent equations 
are constructed as an extension of the non-crossing approximation (NCA) to include the vertex 
part yielding the bound state. A numerical calculation demonstrates the excitonic bound state 
inside the semiconducting gap for single-particle excitations, and a discontinuity at the gap 
edge for magnetic excitations. 

KEYWORDS: exciton, bound state, mixed-valence semiconductor, c-f Coulomb repulsion, single-particle 
excitation, dynamical magnetic susceptibility, non-crossing approximation (NCA) 



1. Introduction 

A class of materials called the Kondo semiconductor or 
the mixed-valence semiconductor has attracted much at- 
tention by its peculiar nature; its semiconducting behav- 
ior is characterized by a small energy gap emerging grad- 
ually with decreasing temperature. It indicates strong 
correlations among electrons as the origin of the energy 
gap. 

Typical compounds classified into this category include 
SmB6, YbBi2 and TmSe. One of noticeable features in 
the semiconducting phase is magnetic excitations de- 
pending strongly on temperature. Neutron scattering ex- 
periments on SmP>6 have shown that an inelastic peak ap- 
pears around 14meV only in the semiconducting phase. 1 
With increasing temperature, the peak disappears cor- 
responding to a transition to the metallic phase. Fur- 
thermore, the intensity of the excitation depends much 
on the wave-vector. The peak with strong temperature 
and wave-number dependences cannot be attributed to 
crystalline electric field excitations. On the other hand, 
experiments on YbBi2 have shown several narrow inelas- 
tic peaks near the spin-gap edge. 2-4 Intensity of one of 
the peaks increases strongly with decreasing temperature 
as in SmB 6 . 

The above-mentioned magnetic excitation have been 
ascribed to an exciton, which here means a bound state 
localized around the rare-earth site. 5-7 The exciton is 
originally a general concept in the optical absorption in 
semiconductors. Due to the long-range Coulomb interac- 
tion, the particle and hole excited by a photon attract 
with each other to form the hydrogenic bound state. As 
a result, the optical absorption takes place at discrete en- 
ergies lower than the energy gap. Based on this concept, 
such a bound state due to the Coulomb interaction be- 
tween the conduction electron and the rare-earth ion has 
been postulated in the mixed-valence systems. Namely, 
a conduction electron produced by mixing with a 4f elec- 
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tron are trapped around the rare-earth site by the effec- 
tive attraction of the resultant 4f hole, and analogously 
a valence hole are trapped by the effective attractive po- 
tential around the 4f electron. Unlike the optical absorp- 
tion, the excitonic bound state in this case can exist at 
the ground state, since the valence of the rare-earth ion 
may be continuously fluctuating through hybridization. 

Although the exciton has been invoked to interprete 
the characteristics of the semiconducting phase in SmBg 
and YbB 12 , most descriptions in mixed- valence semicon- 
ductors are phenomenologically given on the assumption 
of the presence of the excitonic bound state. From the 
microscopic point of view, it is not obvious whether the 
bound state does appear in mixed- valence compounds. 
The objective of this paper is to clarify an existence of 
the excitonic bound state in the mixed-valence semicon- 
ductors on the basis of a microscopic model. 

Since the excitonic bound state is due to a single-site 
effect, it is reasonable to consider a single rare-earth im- 
purity embedded in a semiconductor. We restrict our at- 
tention to the single-impurity model. We consider the 
Coulomb repulsion between the 4f electron and conduc- 
tion electrons as a source of the bound state, in addition 
to interactions in the conventional Anderson model. 

Concerning the low-energy dynamics, it is known that 
an influence of the c-f Coulomb repulsion can be ab- 
sorbed into other parameters. In fact, a numerical renor- 
malization group (NRG) calculation has shown that the 
single-particle excitation spectrum computed with the 
c-f Coulomb term can be obtained from another set of 
parameters without it. 8,9 No qualitative difference has 
arisen by the c-f Coulomb interaction in the metallic 
phase. However, the results cannot be applied simply to 
the semiconducting phase. 

In this model, the Kondo effect occurs provided the 
Kondo temperature exceeds the energy gap. 10 ' 11 On the 
other hand, the Kondo effect is suppressed if the en- 
ergy gap is large compared to the Kondo temperature, 
on which we concentrate our attention. To construct 
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an approximation which is proper in the semiconduct- 
ing phase, we perform a perturbation expansion from 
the atomic limit. We shall take account of leading con- 
tributions for the c-f Coulomb interaction as a vertex 
correction to hybridization. Then, an appearance of the 
excitonic bound state will be demonstrated in the single- 
particle and the magnetic excitations. 

This paper is organized as follows. In the next section, 
we write the Anderson model with the c-f Coulomb repul- 
sion in terms of A-operators, and give a brief summary of 
the resolvent method, by which we perform a perturba- 
tion analysis in the following section. Section 3 discusses 
a vertex part leading to the excitonic bound state. The 
binding energy is examined, in detail, for the on-site in- 
teraction in §4. A self-consistent perturbation theory in- 
cluding the vertex part are developed as an extension of 
the NCA for the model with the c-f Coulomb interaction 
in §5. Numerical results are given for the single-particle 
and the magnetic excitations in §6. We conclude with a 
summary in the final section. 

2. Model and the Resolvent Method 

We start from the Anderson model, which expresses a 
local nature of 4f electrons and strong correlations among 
them. The Anderson model is known to bear character- 
istics of the Kondo effect and the valence fluctuations. 9 
We shall consider, in addition to the conventional Ander- 
son model, the Coulomb repulsion between 4f electrons 
and conduction electrons to capture properties emerging 
from the local bound state. 

We assume the strong Coulomb repulsion between 4f 
electrons, and thus restrict 4f states to 4P and 4P +1 . In 
order to describe such a situation, we write the Hamilto- 
nian in terms of A-operators defined by A 7 y = |7)(7'|, 
where I7) denotes 4f™ and 4P +1 states. On the other 
hand, |a) stands for 4P states, and |/3) 4P +1 states 
hereafter. In this notation, the completeness of 4f states 
within the restricted Hilbert space of 4P and 4P +1 is 
represented as 

yjXy 7 = y]x aa + y^xpp = 1. (1) 

7 a /3 

With use of the A-operators, the Hamiltonian is written 

as 

H = H a + i?h y b + H c f , (2) 

where 

Hq = 22 e k c km Ckm + E 1 X 11 , (3) 

fern 

ffhyb = V {P\fi\a)X 0a c km + h.c. (4) 

kma{3 

Here ek denotes the energy of conduction electrons with 
respect to the chemical potential, and bases of the con- 
duction electrons are represented by the Wannier orbital 
centered at the impurity site. E y is the energy level of 
4f states I7) in the atomic limit, and only their relative 
values are relevant. We take V as a real constant inde- 
pendent of k. 
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Fig. 1. Diagrammatical representations of the c-f Coulomb inter- 
action corresponding to each term of eq. (6). The solid, dashed 
and wavy line represent conduction electron, 4f n+1 state and 4f™ 
state, respectively 

The term H c f in eq. (2) represents the Coulomb repul- 
sion between conduction electrons and 4f electrons, and 
is written as 

H c! = ^ U cl {q)cl +qm c km {N[ - v), (5) 

kqm 

where Nf is the 4f number operator. We have introduced 
a constant z/, which produces a potential scattering inde- 
pendent of 4f electrons, to determine a standard for the 
c-f interaction. For 4f°-4f 1 systems, v is, in most cases, 
taken to be 0. In the general case of n, however, it is 
optional between n and n+ 1. The term of v gives rise to 
a change of the local density of states (DOS) of conduc- 
tion electrons. Hence, we shall choose either n or n + 1 
for the value of v so as to largely cancel out the Hartree 
term. The wave-number dependence of the Coulomb in- 
teraction U (q) is kept for an analytical analysis in §3. In 
constructing self-consistent equations in §5, the interac- 
tion is eventually treated as a point contact as is often 
done in the actual evaluations. In the restricted Hilbert 
space of 4P-4P +1 , H c { is rewritten in terms of the X- 
operators as 

kqm (3 

-J2 u o {q)c{ +qm c km Y . (6) 

kqm a 

where we have introduced (/-dependent coupling constant 
Ui(q) = (n + 1 - v)U cl (q) and U Q (q) = (y - n)U cf (q). 
Figure 1 shows diagrammatical representation of the c-f 
interaction. 

In order to deal with the c-f Coulomb interaction, 
it is convenient to begin with the atomic limit. We 
can utilize conventional approaches developed for the 
Kondo model or the Coqblin-Schricffer model. Resolvent 
method is one of frameworks to perform the perturba- 
tion expansion from the atomic limit. 9, 12-1 A similar 
approach is also possible by using Abrikosov's pseudo- 
particle method. 15-17 In this paper, we make use of the 
resolvent method. 

The partition function Z{ of the 4f part is given in 
terms of the resolvent as 

Zt= I 5se-*£ A rW' (7) 
Jc 27ri 7 

where the contour C encircles all singularities counter- 
clockwise. The resolvent is given by 

Ry(z) = [z - E 1 ~ ^{z)]- 1 . (8) 
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The self-energy part E-y(z) is evaluated by the perturba- 
tion theory. The perturbation term is i?h y b + H c { in our 
model, eq. (2). 

Once the resolvents are obtained, physical quantities 
such as the single-particle excitations and the response 
functions are evaluated with use of them. Since expres- 
sions of physical quantities always involve the Boltzmann 
factor e~* z , the resolvents are dealt with at real frequen- 
cies in actual computations. 

We should note here that the term i?hyb can play a 
role similar to H c f. Namely, the second order process 
with respect to hybridization is regarded as an effec- 
tive exchange or Coulomb-type interaction. The inter- 
action naturally becomes antiferro-exchange or repulsive, 
as in the Kondo and the Coqblin-Schrieffer model derived 
from the Anderson model. We should therefore bear in 
mind that the repulsive c-f interaction may exist even 
though H c f is missing. Such a situation will be demon- 
strated by numerical calculations in §6. 

3. Intermediate State Interaction in Valence 
Fluctuations 

In the optical absorption in semiconductors or insula- 
tors, final-state interaction, that is, the Coulomb attrac- 
tion between excited particle and hole, leads to the bound 
state called exciton, which is located in the energy gap. 
The relative motion of the particle-hole pair is described 
by the hydrogenic wave function with s-symmetry. In the 
Green function method, such bound states are obtained 
by summing ladder-type diagrams in the perturbation 
series for the response function. 18 ' 19 

Similar process is possible in the mixed-valence sys- 
tems, where 4f states are continuously fluctuating be- 
tween 4P and 4P +1 . Unlike the optical absorption, such 
a process takes place between the intermediate states 
of hybridization. Accordingly, the bound state should 
emerge into a variety of dynamical quantities related to 
the 4f electrons. To describe the dynamics, we begin with 
the NCA, which is the self-consistent perturbation the- 
ory with respect to hybridization. 20 ' 21 The NCA gives a 
good description of the feature of the valence fluctuation 
apart from the c-f Coulomb interaction. The non-local 
c-f Coulomb interaction will be taken account of as a 
vertex correction to the NCA self-energy. We take, as 
the vertex part, the intermediate- state interaction of the 
ladder-type diagrams to derive an excitonic bound state, 
following refs. 18 and 19. 

In this section, we omit the spin and orbital degener- 
acy for simplicity, since their degrees of freedom are not 
relevant to the following discussion on formation of the 
excitonic bound states. We assume both 4P and 4f n+1 
states to be singlet. 

3.1 Interaction between 4P +1 state and holes in valence 
band 

We first consider the self-energy part £ n (z) of the 
4f™ state. Figure 2 shows typical diagrams of order V 2 . 
Figure 2(a), which does not include the c-f Coulomb in- 
teraction, yields all diagrams without crossing of conduc- 
tion lines when the resolvent is regarded as the renormal- 
ized one. In addition to the second order process with re- 
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(a) (b) 




Fig. 2. Examples of the self-energy part S n (z) of 4P state of 
order V 2 . 



(a) (b) 




Fig. 3. Sketches of the (a) valence and (b) conduction band con- 
sidered in subsection 3.1 and 3.2, respectively. 

spect to hybridization, diagrams in Figs. 2(b) and (c) in- 
clude interactions between 4P +1 intermediate state and 
the excited electron or hole. The solid lines going in the 
same direction as dashed line represent electron propaga- 
tors, and correspond to the distribution function of the 
hole. On the other hand, lines with the opposite arrow 
stand for holes, and correspond to the electron distribu- 
tion function. We take almost filled valence band located 
E g below the chemical potential as shown in Fig. 3(a). 
In this situation, a contribution from the diagram in 
Fig. 2(c) is negligible because of the hole distribution 
function, provided that temperature is much lower than 
E g . Consequently, the ladder- type diagrams like Fig. 2(b) 
give leading contribution among all intermediate state 
interactions in the semiconducting state. 

The contribution from the ladder-type diagram modify 
the NCA self-energy as follows: 

E n (z) = V 2 f(tk)Rn+i {z + e k )A n+1 {k, z), (9) 

k 

where /(e) is the Fermi distribution function. We have 
introduced the dimensionless vertex function A n+1 (fc, z). 
The vertex part incorporates consecutive scattering of 
holes by 4f 1+1 state, and is given by the following integral 
equation: 

A n+1 (k,z) = 1-J2 UtQk' - fc|)/(e fe ') 

k' 

x R n+l (z + e k i)k n+l (k' ,z). (10) 

Diagrammatical representations of eqs. (9) and (10) are 
shown in Fig. 4(a) and (b), respectively. The minus sign 
of the second term in the right-hand side of eq. (10) 
comes from one permutation process of the conduction 
operators. It indicates the fact that the interaction be- 
tween electron and hole is attractive. When V depends 
on the wave vector k, the following argument is also avail- 
able by defining another vertex function A„ + i(fc, z) with 
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Fig. 4. Diagrammatical representation of (a) the self-energy part 
E n (z) of the 4f™ state and (b) the equation of the vertex part 
A n _|_i(fc, z) corresponding to eq. (10). 



dimension of energy. Then an equation for A n+1 (fc, z) is 
given by replacing the first term on the right-hand side 
of eq. (10) with V h . 

In Fig. 4(a) and (b), we could represent the vertex 
part in a different manner as shown in Fig. 4(c). Both 
diagrams express an identical function since now V is 
independent of k. We shall use both representations to 
construct a self-consistent scheme in §5. 

We solve eq. (10) assuming a simple situation, and 
then confirm that the vertex part introduced above can 
give a description of the excitonic bound state. We con- 
sider only the valence band with a parabolic energy dis- 
persion of the form: 

e k = -E s - k 2 /2m. (11) 

The valence band is assumed to be completely filled, i.e., 
/(e/c) = 1) which is valid at low temperatures such that 
T -C E g is satisfied. Then, R n+ i(z) can be replaced by 
bare one as follows: 

R n+1 (z)^R { °i 1 (z) = {z-A}-\ (12) 

where A = E n+ \ — E n denotes an energy difference be- 
tween 4F 1 " 1 " 1 and 4P states. The above replacement is 
now exact since there is no unoccupied state to mix with 
4P +1 state. Correspondingly, eq. (9) gives the exact self- 
energy. 

To derive the excitonic bound states, it is reasonable to 
deal with the equations in the real space. We introduce 
an auxiliary function S(r, z) defined by 



S(r,z) 



k 



e «*"-J?JW 1 (z + efc )A Tl+1 (fc,z) 
A n+ i(k,z) 



— A — E e 



fc 2 /2m' 



(13) 



which is the Fourier transform of the integrand of eq. (9) . 
The self-energy part £„(z) is represented with use of 
S(r, z) as follows: 

X n (z) = V 2 S(0,z). (14) 
Operating (z — A — E s + V 2 /2m) on eq. (13), we obtain 
[z - A - E e + V 2 /2m]5(r, z) = ^ e ik - r A n+1 (k, z) 



6(r) 



C/iM5(r,z), 
(15) 



where we have used eq. (10) in deriving the final expres- 
sion. This equation includes the one-body Hamiltonian 
i? cx = — V 2 /2to — Ui(r) with attractive potential, which 
can yield a bound state depending on U\(r). The bind- 
ing energy will be examined in detail for the short-range 
interaction in the next section. We represent eigenvalues 
of H cx by ej as follows: 

H ex 4>i(r) = [-V 2 /2m - E7i(r)]&(r) = e^r). (16) 

We then expand S(r, z) in terms of the orthonormal 
bases (f>i(r). Solving eq. (15) with respect to S(r,z), we 
obtain 

v mm?) 



S(r,z) 



(17) 



This expression is substituted for eq. (14) to yield 

l<M0)l 2 



Z n (z) = V 2 J2 



(E g + ei) 



(18) 



This expression describes a virtual excitation of a hole, 
whose energy is modified by the intermediate-state in- 
teraction into E g + ei. The energy ej given by eq. (16) 
may be negative due to the effective attraction, describ- 
ing a bound state localized around the 4f site. Note that, 
because of the factor |</>i(0)| 2 , only s-wave parts among 
4>i{r) contribute to the self-energy. This results from the 
constant V, where a spread of the 4f electron is neglected. 

3.2 Interaction between 4P state and electrons in con- 
duction band 

The self-energy part S n+ i(z) of 4P +1 state can be 
constructed in a manner similar to E„(z). In this case, 
interactions between 4f intermediate state and excited 
conduction electrons are taken into account. In the ladder 
approximation, the NCA self-energy is altered into 

En+iC*) = V 2 f(- e k)Rn(z - e fe )A„(fc, z), (19) 

k 

where the vertex function A„(fc, z) is defined by 
A n (k,z) = l-J2M\k-k'\)f(-e kl ) 



x R n (z - e k ')A n (k', z). 



(20) 



As in the A„ + i(fe, z), the minus sign of the second term 
of the right-hand side indicates an attractive interaction. 
Corresponding diagrams are shown in Fig. 5. The vertex 
part A n (fe, z) as well as A n+ i(fe, z) can be expressed in a 
couple of ways as shown in the Fig. 5(c). 

We solve eq. (20) for a toy model following the previous 
analysis. We suppose a conduction band with a parabolic 
dispersion represented by 



e k = Eg + fc 2 /2m, 



(21) 



which is sketched in Fig. 3(b). We assume T <C E s so that 
the conduction band is taken to be completely empty, i.e., 
f(~ e k) — 1- Then, R n (z) can be replaced by bare one 

R n (z)^RW(z)=z-\ (22) 
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(c) 
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Fig. 5. Diagrammatical representation of (a) the self-energy part 
>;„ . i (z) of the 4f™+ state and (b) the equation of the vertex 
part A n (fe, z) corresponding to eq. (20). 



without any approximations. At the same time, the self- 
energy part S„+i(z) given by eq. (19) becomes exact. On 
the above assumption, eq. (20) is solvable in the same 
procedure as the derivation of eq. (18). We show only 
the result for £„ + i(z): 



^n+l{z)=V 2 



^ 2 5>K0)| 2 i?£ 0) (^-£ g - e ;:), (23) 

i 

where ^(r) and are determined by 

KA(r) = [-V 2 /2m - U (rM(r) = e'^r). (24) 

Equation (23) indicates a virtual excitation of the elec- 
tron with energy E s + e[ , which may be smaller than E g 
due to the effective attraction. Accordingly, it describes a 
situation that the 4f electron excited into the conduction 
band is bound around the impurity site. 

4. Binding Energy of the Local Bound State 

It has been shown that, in the absence of the ther- 
mal population of electrons or holes over the energy gap, 
virtual excitation processes are governed by the one- 
body Schrodinger equation with the attractive potential, 
eqs. (16) and (24). It is apparent that the hydrogenic 
bound states appear provided the interaction extends to 
long range. In the metallic compounds, however, it is not 
clear because the interaction is short-ranged due to the 
strong screening effect. In order to see the condition for 
occurrence of the bound state, we consider an attractive 
<5-function potential for each dimension 

-Ua d S d (r), 



U{r) 



(25) 



where a denotes the lattice constant and d is the spacial 
dimension. 

From the Schrodinger equation with the attractive po- 
tential, we obtain the following condition which deter- 
mines the binding energy E^: 



(26) 



where k is defined by k — \/—2mEh and N denotes a 
number of sites. Since we are now interested only in the 



bound state, we have restricted E\> to negative. We intro- 
duce a cutoff energy D and corresponding wave number 
k c = \/2mD so as to satisfy the relation 



/ dep(e) = p D. 
Jo 

Then p(e) is given by 



dp ( e \ d/2-i 



(27) 



(28) 



for each dimension d. 

For one dimension d 
eq. (28) to yield 

PoUk c 
1 = arctan 



1, eq. (26) is evaluated with 
np Uk c 



2k 



(29) 



where we have used the approximation arctan(fc c //i) 
7r/2. Hence Eh is given by 



E h = — 



n 2 D 



(p Uf 



(30) 



We note here a resemblance to the bound state appear- 
ing in the energy gap of superconductors by magnetic 
impurities. 22 Since the DOS of the quasi-particles in su- 
perconductors has e _1 / 2 -dependence at the gap edge, the 
binding energy is given by a formula similar to eq. (30). 
For two dimensions d — 2, eq. (26) leads to 

E b = -D(e 1/poU - I)' 1 . (31) 

When e 1 ' /paU ^ 1 is satisfied, which is a plausible condi- 
tion in real systems, eq. (31) is reduced to 

E h = -De- 1/poU . (32) 

This expression is identical to the well known formulae 
for the Kondo temperature Tk and the superconducting 
transition temperature. It seems that the expression like 
eq. (32) is a common form emerging from the constant 
DOS. 

In the case of three dimensions d = 3, in contrast to 
the lower dimensions, there exists a threshold value of U. 
A bound state appears only when U is large enough so 
that poU > 1/3 is satisfied. Then Eh is given by 

E, = -~ ( i - ~ ~77 } , (33) 



1 

3put7 
applied the 



approximation 



where we have again 
arctan(/c c /«;) ~ n/2. 

It turns out, from the above discussion, that there 
may exist one excitonic bound state even though the c- 
f Coulomb repulsion is on-site. As expected, the bound 
state becomes more stable, the larger the DOS of con- 
duction electrons exists near the gap edge. In the semi- 
conducting state, the energy gap weakens the screening 
effect because of an absence of the particle-hole excita- 
tion with an infinitesimal energy transfer. Accordingly, 
the c-f interaction is not necessarily on-site there. Hence, 
there is a possibility that plural bound states exist in 
semiconducting phase. 
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(a) 



Fig. 6. Diagrams representing the generating functional <E> which 
includes the consecutive scattering due to the c-f Coulomb inter- 
action. The first term represents the NCA. 



5. Self-Consistent Perturbation Theory 

In constructing an approximate scheme by a partial 
resummation of perturbation series, it is important to 
let the scheme preserve the natural characteristics such 
as the analyticity and the sum-rule. It has been shown 
for the system of electron gas that approximate schemes 
constructed from generating functionals $ of the skele- 
ton diagrams conserve the proper characteristics. 23 Such 
an approach is also available for the resolvnent method, 
which deals with the restricted Hilbert space. 20 In this 
section, we derive equations for the resolvent and some 
dynamical quantities within the conserving approxima- 
tion, following the analysis of the vertex part in §3. The 
c-f Coulomb interaction U c f is considered to be indepen- 
dent of the wave number, since the ^-function potential 
has turned out to give a bound state in the previous sec- 
tion. 

5.1 Self- energy part 

The self-energy part is obtained by functional differ- 
entiation of <E> with respect to R~ f (z). Diagram in Fig. 6 
shows the generating functional which leads to the self- 
energy parts including the ladder-type intermediate-state 
interactions, Figs. 4 and 5. The first term generates the 
NCA self-energy, and others modify it so as to incorpo- 
rate terms leading to the excitonic bound state. 

Generating functional in Fig. 6 are reminiscent of 
those introduced in extensions for the finite-J7 Anderson 
model, where the vertex correction incorporates vertual 
excitations to 4? states. 24-26 We can observe a structural 
correspondence between both generating functionals. Ac- 
cordingly the knowledge of the schemes concerning the 
finite-C/ Anderson model are applicable in this model. 

The self-energy parts generated from $ in Fig. 6 are 
diagrammatically represented in Fig. 7. The first term 
is the ladder-type diagram introduced in Figs. 4 and 5. 
In addition, the second diagram is necessary for the con- 
serving approximation. The third diagram subtracts a 
double counting part. We introduce auxiliary functions 
P a {z) and Pp(z) defined by 



P a {z)= / dep(e)f(-e)R a (z-e), 



Pp{z)= / dep(e)f(e)Rp(z + e). 



(34a) 



(34b) 




(b) 




+ 



+ 






Fig. 7. Diagrammatical representation of the self-energy part (a) 
T, a (z) and (b) Sfl(z) in the conserving approximation. 



Since now the c-f Coulomb interaction is on-site, the ver- 
tex function A n _|_i(fc, z) and A„(fc, z) introduced respec- 
tively in eqs. (10) and (20) become independent of the 
wave vector k. The integral equations are consequently 
solvable and lead to 



A a (z) = [l + U Q P a {z)}-\ 
A f3 (z) = ll + U 1 Pp(z)}-\ 



(35a) 
(35b) 



It is clear that the vertex part A 7 (z) reduces to unity 
if the c-f Coulomb interaction does not exist. The self- 
energy part E Q (z) of 4P state |a) is given with use of 
Pp(z) and kp(z) by 

E a (*) = £ WaJp, CO [A/J (*)-!] 



+ / dep(e)f(e)R 0a (z + e) , (36) 



(37) 



Rp a {z) =R(3(z)A 2 a (z), 



where W af} is defined by W a/3 = V 2 £ m |(/3|/tja)| 2 . In 
an analogous way, the self-energy part E^ of the 4P +1 
state |/3) is given by 



(z)[Aa(*)-l] 



dtp(e)f(-t)R a p(z 



R a p{z) = R a (z)A}(z). 



(38) 
(39) 



These equations are solved numerically combined with 
eq. (8). 

We define spectral function ry 7 (e) of the resolvent and 
that of the defect propagator £ 7 (e) . They are expressed 
in terms of the resolvent as 27 



??7 ( £ ) 

£r(<0 



-7r -1 Imi? 7 (e + ), 
Z^e-^i-Tr-^ImR^ 



(40) 
(41) 



where e + denotes e plus an infinitesimal complex number. 
ry 7 (e) stands for intensity of an addition of the state I7), 
while £ 7 (e) a removal of the state I7). Although £ 7 (e) 
analytically relates with the resolvent in eq. (41), the 
Boltzmann factor makes it difficult to compute £ 7 (e) at 
low temperatures. To avoid the difficulty, we transform 
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Fig. 8. Diagrammatical representation of the 4f electron Green 
function. 

the integral equations for the resolvent into a form con- 
venient for the numerical calculation of £ 7 (e) (see Ap- 
pendix A). With use of these spectral functions, dynam- 
ical quantities are represented without the Boltzmann 
factor. £ 7 (e) is further used to evaluate the partition func- 
tion Zf in eq. (7). Thermodynamic quantities such as 
the entropy and the specific heat are evaluated from T- 
derivatives of the spectral functions, which can be com- 
puted with use of another set of integral equations. 28 

5. 2 Single-particle Green function 

We now derive the single-particle Green function from 
the generating functional $ in Fig. 6. The 4f Green func- 
tion is defined by 

Gf, m (ie n ) = - f dr{ff n {r)fl)e^ 
Jo 

a/3 

where e„ = (2n + 1)ttT is the fcrmion Matsubara fre- 
quency and the superscript H indicates the Heisenberg 
picture. G a p(it n ) is defined in terms of the X-operators 
by 

G Q/3 (ie„) = - / dT(X^(T)Xp a )e^\ (43) 
Jo 

We consider a response to an infinitesimal external 
field which changes the 4f number. As far as hybridiza- 
tion process is concerned, our approximation is equiva- 
lent to the NCA, since the vertex parts merely describe 
interactions between the intermediate states. Therefore, 
the formula of the 4f Green function in the NCA is ap- 
plicable in this case. Hence G Q/ g(ie„) is given by 

G aP (ie n ) = i- / p,e- f3z R a (z)R p (z + ie n ), (44) 
Zf J c Ziri 

which is expressed by diagram in Fig. 8. Performing ana- 
lytic continuation ie„ = uj + to real frequencies, we obtain 
the DOS p a f}(uj) = — 7r~ 1 ImG Q , ( g(a; + ) of 4f electrons 

Pa/sH= /de[e a (e)»?/j(e + w) + »7 a (e)^(e + w)]. (45) 

5. 3 Two-particle correlation junctions 

We proceed to the two-particle correlation function. In 
general, 4f electrons possess the orbital degrees of free- 
dom. Correspondingly there are various kinds of corre- 
lation functions other than spin-spin and charge-charge 
ones. Therefore it is reasonable to consider formula for 



z, y Y" 




z+iv„, y' y" 



z, y y'" z, y z, y z+e, V y'" 




z',Y y" z'.y' z',Y z'+s,X y" 



Fig. 9. Diagram of the tow-particle correlation function. Double 
line means all 4f states 7 including 4f n and 41™+ , namely, it can 
be thought of as dashed or wavy line. 

correlation functions in a general form with four sub- 
scripts. At the end of this subsection, we shall apply it 
to the dynamical magnetic susceptibility in a model for 
Ce compounds. 

5. 3. 1 General formula 

Let us consider a correlation function of operators A 
and B which do not change the 4f number, i.e., (a|A|/3) = 
(a\B\(3) — 0. The 4f response function is defined by 

XabK) = / dr{T T A n (r)B)e iv ^ 
Jo 

= E (7l^l7 / )(7"l-B| 7 '")x 7 y )7 "y»K), 

77'7"7"' 

(46) 

where v n = 2mrT denotes the boson Matsubara fre- 
quency and X77',7"7'" (i^n) 1S defined in terms of the X- 
operators by 

X77',7"7"'CO= / dT(T T X^,(T)X Y , y „)e^. (47) 
Jo 

The NCA requires the vertex correction in the two- 
particle correlation function if, for example, both 4P 
and 4f rt+1 states have the magnetic moment 29 or the hy- 
bridization is anisotropic. 30 In the present model, an ap- 
proximation consistent with $ in Fig. 6 needs the vertex 
correction. The 4f response function is given by 

(' \ f ~P Z Y\ ( ' \ 

Zf J c Ziri 

(48) 

n 77 / j7 » 7 »'(z, z 1 ) satisfies the following Bethe-Salpeter- 
type integral equation: 

-L-L'-yy ' ^lyffsyt/t ( Z j Z ^ Z ^ -f?»y I ^ Z ^ ^ Sry-y/f/ S *y ' ry ft 

+ E / der 77',AA' (z, z'; e)n A 'A,7"7'" (z + e, z' + e) > , 

(49) 

where sum of A and A' are taken for all 4f states. 
r 77 / j7 " 7 "/ (z, z'; e) denotes an irreducible vertex function. 
Diagrammatical representations of eqs. (48) and (49) are 
shown in Fig. 9. 
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We classify II 77 <. 7 " 7 »< by their subscripts, i.e., we make 
use of subscript a (4P states) and f3 (4f™ +1 states) in- 
stead of 7 (whole states). It is sufficient to introduce four 
classes of functions: U { ^) aa ,, U^)^,, ^%)^, fj ,„, and 



n^jj a ii a m- Analogously, the irreducible vertex functions 
r 77 ' i7 » 7 '» are classified into four parts: r^j? aa ,, r^?^,, 

pp> 3" 8"" an< ^ W aV" ^ n * ne P resen t approximation, 
each function is given by 

x [A f3 (z ~ e)kp{z' - e) + A a ,{z)A a (z') - 1], (50a) 

x [A /3 /(z)A /3 (z') + A a (z + e)A a ,(z' + e) - 1], (50b) 

= | de'p(e')/(-e'Me + e')/(e + e') 



+ W a/3/ R al3 ,(z' -e')A^(z-e') 



(50c) 



1 aa'a'a\ z i z ' b ) 



-U J de'p(e')f(e')p(e + e')f(-e-e') 
x ^\W a pRp a {z + e')A a < [z + e J 

W Q ^ J R / 3 Q -(z' + e')A Q (z + e')], 



(50d) 



where we have introduced VFa Q ',/3/3' = 
V 2 E m (P'\fUa)(a'\f m \P)- The above ' quantities 
are diagrammatically expressed in Fig. 10(a). The 
first and second irreducible vertices are what should be 
considered in the NCA, namely, they do not become zero 
even if U c i = 0. On the other hand, the third and fourth 
ones are entirely due to the c-f Coulomb interaction, 
that is, they dissappear if U c i = 0. In deriving eqs. (50c) 
and (50d), we have introduced effective c-f Coulomb 
interactions Up(z) and U a (z) corresponding to diagrams 
in Fig. 10(b). Since the c-f Coulomb interaction is now 
on-site, the integral equations are solvable and give 
Up(z) = U x Ap{z) and U a {z) = -U A a (z). 

After analytic continuation \v n = lo + to the real axis 
in eq. (48), we obtain an expression for the magnetic 
spectrum 



Lnx 77 ',yv"(w + ) = (1 - e^") 



1 f de 



Zf J 2vr 



x Re[II 77 / 7 » 7 «/(e , e + u + ) — H 77 / 7 » 7 ///(e + , e + oj + )]. 

(51) 



The static susceptibility x(0) is given by 

X77',7"7"' (0) = "77" I ® iDtJI™/ ~//~/ 
J 7T 



(52) 




V 6 



A E 



z,P 




z+£,p 


— >~ 


p(H) 


r->— 


— <- 

«'.P' 




■-<— 

z'+e, P' 


z, ot 




z-e, a 




p(00) 




z', a 1 




z'-e, a' 




Z'+E', P 



(b) 



--<--( )-<- 

Z+E Z+E' 




Z-E Z-E' Z-E" 



Fig. 10. Diagrammatical representations of (a) the irreducible 
vertex parts vQ?) ,, V ; „„,, T„}) ,,,,, and r' 00 ' , , and (b) 

pp ,ckck ' aa ,pp pp ,p p aa ,a «' v ' 

equations for the effective c-f Coulomb interactions Up{z) and 
U a {z). 



5.3.2 Dynamical magnetic susceptibility in a model for 
Ce compounds 
For the systems of 4f°-4f 1 , only II^ 11 ) contributes to 
the dynamical magnetic susceptibility since 4f° has no 
magnetic moment. In the absence of magnetic field, con- 
tributions from terms including r' 10 ) and F^ 01 > as irre- 
ducible vertices become zero, since sum of the magnetic 
moment appears in such terms. Hence, T^ 11 ) is the only 
irreducible vertex relevant to the magnetic susceptibility. 
Dynamical magnetic susceptibility XM^n) is then given 

by 

> N\C f dz _fl ZTT , , 

XM{Wn) = 7T~ I 7T^ e Tl M (z,z + iiy n ), (53) 

6i Jc 27Tl 

where C is the Curie constant and A^i is the number of 
degeneracy of 4f* states. Hm(z, z') satisfies the following 
integral equation: 

U u {z,z') = R 1 (z)R 1 {z')!l-V 2 U 1 fdep(e)f(-e) 



x [R 01 (z - e)Ai(z' - e) + R i(z' - e)A 1 (z - e)} 

x J de'p(e')f(e')n M (z - e + e' , z' - e + e')|, (54) 
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which follows from eqs. (49) and (50c). Imaginary part 
of Xm(w) and the static magnetic susceptibility xm(0) 
are evaluated by combined with eqs. (51) and (52), re- 
spectively. In computing these equations including the 
Boltzmann factor, we transform eq. (54) into suitable 
forms (see Appendix B). 

6. Numerical Results 

In this section, we show numerical results of the self- 
consistent equations developed in the preceding section. 
We consider a model for Ce compounds, where 4f state 
fluctuates between 4f° and At 1 . Degeneracy of 4f* states 
is assumed to be six (j = 5/2), and we do not consider 
the crystalline electric field splitting. In this model, all 
matrix elements of f m are unity. We take a constant DOS 
of conduction electrons with energy gap 



P(e) = ^{D + E g 



|e|)0(|e|-25 g ). (55) 



We take D as a unit of energy i.e., D = 1. 

Numerical calculations are first performed for param- 
eters where 4f electron is almost localized, to which the 
Ce ions usually correspond. The case where 4f states are 
almost unoccupied will be examined as well for compar- 
ison. 

6.1 4f -dominant case 

As a preliminary to systems with energy gap, we show 
results for the metallic phase to examine the influence 
of U c { on the Kondo effect. Figure 11 shows tempera- 
ture dependences of the entropy Sf and the 4f number 
nf for several values of U c { with ef = E\ — Eq = —0.2 
and V = 0.1. The entropy has been computed without 
numerical differentiation by solving integral equations 
for T-derivatives of the resolvent. 28 It turns out, from 
Fig. 11(a), that the repulsive values of U c { increase the 
Kondo temperature. The fact means that the effect of 
U c { is, at least in the metallic phase, simply regarded 
as a increase of V. Correspondingly, the 4f number de- 
creases with increasing U c { as shown in Fig. 11(b). The 
DOS Pf(w) of 4f electrons near the Fermi level is shown 
in Fig. 12. The Kondo resonance shifts toward the higher 
energy indicating increase of the Kondo temperature. 
These results are consistent with the previous study by 
means of the NRG. 8 

We now proceed to the case of semiconducting phase. 
Figure 13 shows the DOS p f (w) = -Tr^ImGf (u>+) of 4f 
electrons for E g — 0.04. We vary U c f with other param- 
eters fixed at e f = -0.2, V = 0.1, v = 1 and T = 0.01. 
In these parameter sets, the Kondo temperature is esti- 
mated to be much lower than the energy gap. We can 
observe two peaks of completely different features: broad 
one at about ef and considerably sharp one of about T 
in width being inside the energy gap. The latter is an 
excitation from the excitonic bound state discussed in 
§3. The bound state becomes stable with increasing U c { 
as expected. Correspondingly, a weight of the peak in- 
creases. We here comment upon the parameter v. In the 
case where v is set to be zero with other parameters fixed, 
the bound state shifts to upper energy. Since now the 4f 
number is more than 0.9, we recognize that the Hartree 
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Fig. 11. Temperature dependences of (a) the entropy Sf and (b) 
the 4f number nf for several values of U c t (u = 1) without the 
energy gap in the conduction band. Other parameters are ef = 
-0.2 and V = 0.1. 



term enhances the effect of U c t to yield a larger binding 
energy. 

It is confirmed from Fig. 13 that the exciton state ap- 
pears even though i/ cf — 0. This is because the second or- 
der process with respect to hybridization effectively acts 
as a substitution for the repulsion between 4f and con- 
duction electrons, as mentioned in §2. Furthermore, the 
large degeneracy Ni = 6 of 4f states enhances the effec- 
tive repulsion. Namely, the almost localized model yields 
7Viy 2 /|ef| as a strength of the effective c-f interaction. 
In the corresponding parameters, it is estimated to be 
N\V 2 /\ef\ = 0.3, and therefore leads to the binding en- 
ergy E h ~ —1.3 x 10~ 3 from eq. (32). The estimation 
turns out to be in good agreement with the binding en- 
ergy measured in Fig. 13. 

In order to see the origin of the bound state, we ex- 
amine the spectral functions ?? 7 (e) of the resolvents and 
the defect propagators £ 7 (e). Figure 14 shows the results 
with Ud = 0.3. In general, ?7 7 (e) is finite only above a 
threshold energy at zero temperature, while £ 7 (e) is finite 
below the threshold. We have adjusted an origin of the 
energy for the spectral functions so as to set the thresh- 
old at e ~ 0. It turns out from Fig. 14(a) that the 4i x 
states are located at lower energy than the 4f° state. The 
fact indicates that the ground state is magnetic, and ex- 
cludes the possibility of the Kondo effect as an origin of 
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Fig. 12. The DOS Pf(u>) of 4f electrons near the Fermi level. Tem- 
peratures are much lower than Tk and other parameters are same 
as in Fig. 11. 
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Fig. 13. The DOS pf(w) for several values of U c { with v = 1. 
Other parameters are e { = -0.2, V = 0.1, Eg = 0.04 and T = 
0.01. 



the sharp peak in pf(w). This is also confirmed by the 
residual entropy of Sf = In 6 and the Curie law behavior 
of the static magnetic susceptibility (no figure). Com- 
parison of 7/ 7 (e) with /Of (w) in Fig. 13 indicates that the 
structure of the single-particle excitation is mostly de- 
termined by 770 (e) . It follows that the bound state comes 
from the 4f° resolvent, where the intermediate 4f x states 
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£/£> 




-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 
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Fig. 14. (a) Spectral functions r/ 7 (e) of the resolvents and (b) 
those of the defect propagators £(e) for U c { = 0.3 with the same 
parameters as in Fig. 13. An origin of the energy has been set to 
an arbitrary value (see text). 



interact with the holes in valence band by the effective 
attraction. The bound state appearing in 770(e) remains 
in £o( e ) due to the thermal population. 

An influence of the impurity on the electrons of the 
conduction and valence band appears in the site-diagonal 
part of the Green function. The local DOS p c {e) at the 
4f site is given by 



p c (u)) = Im[g(w~ 

7T 



)+g 2 (u J +)V 2 G l (u J +)l 



(56) 



where g(uo) is the site-diagonal part of the unper- 
turbed Green function for conduction electrons defined 
by g(oj + ) = Y.k9k{^ + )/N. Figure 15 shows p c (w) for 
the same parameters as in Fig. 13. We can see an obvi- 
ous difference from the bare one p{u). Sharp peak is per- 
ceived at the same energy as the bound state observed 
in pf(cu) in Fig. 13. The state inside the gap is caused by 
the impurity scattering represented by the second term 
in cq. (56). 

We next examine an influence of the bound state on 
the magnetic excitations. Figure 16 shows ImxM(w) f° r 
the parameters identical with above. The vertex correc- 
tion for Xm(w) is not necessary to be calculated since now 
v = 1 (then Ui =0). There exist mainly the following 
three structures. One is the quasi-elastic scattering with 
extremely narrow width. This peak is due to the degen- 
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Fig. 15. The DOS p c (u>) of the site-diagonal part of the conduc- 
tion electrons with the same parameters as in Fig. 13. 




0.04 0.06 
(O/D 

Fig. 16. The magnetic spectrum Imxivi(<^) for the parameters in 
Fig. 13. 



eracy of 41" 1 states. In addition, two structures are ob- 
served for inelastic excitations. The peaks are due to the 
edges of the gap observed in pt{u) in Fig. 13. The excita- 
tion energy, or the threshold energy of the excitation, are 
represented as E g ± Wb, where iv^ < denotes an energy 
where the exciton peak is located in pf (ut). is related 
to the binding energy Ey, as E^ = — E g + \co^\ < for the 
4f localized case. The upper inelastic peak at E g + |ajh| 
corresponds to an excitation from the bound state to un- 
occupied 4f states. On the other hand, the lower one at 
E g — is due to a transition between the bound state 
and the occupied states. 

As demonstrated above, the bound state clearly ap- 
pears in the dynamical responses at low temperatures. 
The next interest is their thermal evolution. Figure 17(a) 
shows temperature dependence of pt{un) for U c f = 0.3. 
As temperature increases, the peak of the bound state 
broadens, and eventurally merges with the broad hill 
at T = 0.025, which is close to the binding energy. At 
T = 0.100 a trace of the bound state is seen at the gap 
edge of the combined spectrum. At higher temperature, 
the spectrum is almost composed of the excitation from 
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Fig. 17. Temperature dependences of (a) pf(u)) and (b) Imxivi(<^) 
for U c { = 0.3. Other parameters are same as in Fig. 13. 
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Fig. 18. Temperature dependences of the 4f number ri{ for the 
parameters in Fig. 13. 



the original 4f 1 states. Corresponding to the evolution 
into broadened spectra with increasing temperature, the 
4f number rif decreases as shown in Fig. 18. At low tem- 
peratures, U c { reduces rif in a manner similar to the gap- 
less case. 

We next show temperature dependences of Im^M (w) in 
Fig. 17(b). As for the quasi-elastic peak, it monotonously 
increases and narrows with decreasing temperature. Be- 
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Fig. 19. The DOS Pf(ui) for several values of U c f and v. Other 
parameters are e f = 0.06, V = 0.1, E e = 0.04 and T = 0.005. 



tween the two inelastic peaks at T = 0.010, the lower 
one is reduced at lower temperature T = 0.006, while 
the upper one keeps the intensity. It is conceivable that 
the lower excitation is enhanced due to the thermal pop- 
ulation at temperatures comparable with the energy gap. 

6.2 4f '-dominant case 

We move on to the parameters where 4f states arc 
mostly unoccupied. This situation is in fact far from ac- 
tual Ce ions, but is taken for comparison with the prac- 
tical model. Figure 19 shows the DOS Pf(u>) for e{ = 0.06 
and T = 0.005. Several pairs of values are chosen as pa- 
rameters U c { and v. In the absence of the term of U c {, 
no distinct bound state can be found in contrast to the 
4f localized case. Since the degeneracy of 4f states does 
not work in the second order process of hybridization in 
this case, the effective c-f interaction remains small. For 
the finite values of U c f, we have chosen both v = and 
v = 1 for comparison. In Fig. 19, finite values of U c f with 
v = slightly enhance the binding energy, while larger 
increases are seen in the case of v — 1 . Since 4f number 
rif is about 0.15 to 0.20, the Hartree term is largely can- 
celed with v = 0, but not with v = 1. It follows that the 
Hartree term contributes to the effective c-f interaction 
as in the localized 4f model. 

Figure 20 shows Imv M (w) for the same parameters. 
The vertex correction has been computed for v = 0, 




Fig. 20. ImxM(<*>) for the same parameters as in Fig. 19. 



which gives a finite value of U\. Unlike the previous case, 
no quasi-elastic scattering is found since the ground state 
is a singlet. Inelastic scattering occurs above 2E g = 0.08, 
which correspond to a transition across the energy gap. 
We interpret from the knowledge of pt{u>) in Fig. 19 that 
the spectrum reflects a distribution of the unoccupied 
4f states shown around ui ~ 0.1 in pf(av). By comparing 
spectra for U c { = and U c { = 0.3 [y = 1), we con- 
clude that the bound state influences the threshold of 
the excitation gap. The spectrum for U c f = 0.3 (u = 1) 
is discontinuous at the threshold, which shifts to lower 
energy due to the bound state, while for U c f = the 
inelastic spectrum continuously arises at 2E g . 

7. Summary and Discussion 

We have demonstrated an appearance of the excitonic 
bound state in the semiconducting phase based on the 
extended Anderson model with the c-f Coulomb inter- 
action. The bound state is ascribed to an effective at- 
traction between the 4f electron and valence holes, and 
between the 4f hole and conduction electrons. An im- 
portance of these processes is clarified by the pertur- 
bation theory which incorporates leading contributions 
of the c-f Coulomb interaction as a vertex correction to 
hybridization. Our microscopic derivation of the bound 
state supports previous discussions on the exciton in the 
mixed- valence semiconductors. 1, 5 ' 6 

Compared to the exciton in ordinary semiconductors, 
an important difference is that the bound state can exist 
at the ground state because the 4f electron mixes with 
conduction electrons without optical absorption. In other 
words, a source of the excitonic bound states appears as 
an effective potential in the single-particle Green func- 
tion, while in the case of the optical absorption a direct 
attraction in the two-particle response function leads to 
the bound state. Consequently, the excitonic bound state 
shows up in the 4f DOS as an extra 5-function-like peak 
located inside the energy gap, in addition to a broad 
peak corresponding to the original 4f level. On the other 
hand, if the original 4f level is close to the Fermi level, 
the spectrum becomes single peak structure inside the 
energy gap. 
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In ordinary semiconductors with electron donors or ac- 
ceptors, the attractive Coulomb potential around the im- 
purities gives rise to the bound state. The donor level lies 
below the conduction band and provides electron carri- 
ers by thermal excitations, while the acceptor level lies 
above the valence band and provides hole carriers. On 
the other hand, in the present model with rare-earth im- 
purities, the bound state located around the 4f site do 
not contribute to the electrical conduction. Namely, the 
effective c-f attraction simply produces the 4f level inside 
the energy gap with the chemical potential fixed, on the 
contrary to the donor or acceptor level. 

A source of the bound state in mixed-valence systems 
is not only the Coulomb interaction but also the effective 
interaction due to the higher order process of hybridiza- 
tion. In fact, it has been demonstrated for the well- 
localized 4f electron systems that the excitonic bound 
state appears in the 4f DOS even though U c f = 0. There- 
fore, the bound state naturally appears around Ce ions 
in the well-localized 4f systems as long as an energy gap 
exists. In the case where the 4f states are almost un- 
occupied, however, a distinctive peak by the excitonic 
bound state does not emerge without U c f. These con- 
trastive results are well understood by considering each 
hybridization process as follows. When the 4f electron is 
removed from the localized 4f electron system, the exci- 
tonic bound state take places in the resultant 4f° state, 
whose self-energy is enhanced due to the A^-fold degen- 
eracy of 4f x states. In the opposite situation, the degen- 
eracy never works in the 4f 1 -state self-energy. 

It is interesting that the energy of the excitonic bound 
state |-Et>| and the Kondo temperature Tk are given by 
analogous formulae. Suppose 4f electron is well localized 
and U c f = 0, both the effective c-f repulsion and the 
antiferro-exchange coupling are given by NiV 2 /\ef |, and 
therefore \Eb\ and Tk become identical. When the sys- 
tem is metallic, the 4f electron forms the Kondo singlet, 
whose energy gain is given by Tk determined from the 
above coupling constant. This is the same situation in 
semiconductor, provided the energy gap is much smaller 
than Tk . On the other hand, if there exists a large energy 
gap compared to Tk, the excitonic bound state arises at 
|£b| — Tk inside from the gap edge in pf(oj), instead of 
the Kondo resonance peak. 

It has been shown that the bound state affects mag- 
netic excitations. The situations are different depending 
on whether the gound state has the magnetic moment 
or not. When the system has no magnetic moment at 
low temperatures, low-energy excitations does not take 
place. If there is no excitonic bound state, excitations 
gradually start to occur at the gap energy. In the case of 
the presence of the bound state, on the other hand, the 
threshold of the magnetic excitation becomes discontinu- 
ous. An energy of the threshold is lower than the original 
gap due to the existence of the bound state. For the mag- 
netic ground state, the quasi-elastic scattering appears, 
and as in the singlet ground state the excitonic bound 
state make the spectrum discontinuous at the threshold. 

We have devoted ourselves to the single-impurity 
model in this paper. The system corresponds to such 
a situation where magnetic impurities like Ce are doped 
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into semiconductors. On the other hand, present results 
can contribute to study on the exciton in the periodic 
systems as well. Suppose the excitonic bound state is 
formed at each site, it transfers from site to site through 
conduction electrons. Propagation to other site may lead 
to, for instance, broadening of the local bound states and 
wave-number dependence of the intensity and the bind- 
ing energy. The situation becomes more complicated in 
the Kondo semiconductors, where 4f electrons themselves 
take part in a formation of the energy gap. It is an in- 
teresting issue whether 4f electrons can simultaneously 
play roles both in an evolution of the hybridization gap 
and in the appearance of the excitonic bound state. The 
periodic systems can be dealt with by the present theory 
combined with the dynamical mean field theory. They 
are left for future studies. 
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Appendix A: Equations for the Defect Propaga- 
tor 

The spectral function of the defect propagator £ 7 (e) 
relates with the resolvent analytically in eq. (41). How- 
ever, the Boltzmann factor makes it difficult to per- 
form direct computations at low temperatures. In order 
to obtain them by numerical calculations, we transform 
the whole integral equations for the resolvent, follow- 
ing ref. 27. We introduce an operator defined by 
V^R(e) = -Z f T 1 e-P z iT- 1 ImR(e+). The operator 7>«) 
converts eq. (8) into 

£ 7 (e) = |i? 7 (e+)|V 7 (e), (AT) 

where tr 7 (e) denotes cr 7 (e) = P^E 7 (e). Operating 
on eqs. (36) and (38), we obtain 

«7 a (e) = ]T W aP \p;{e)\ fj (e) + 7r /3 (e)[A /3 (e) - 1] 

+ J deV(e')/(-e')l/ 3 a(e + e')}, (A-2) 
= Y. W ^{ P «(.£)K{e) + n a (e)[A a (e) - 1] 

+ J de'p(e')f(e^ aP (e-e')y (A-3) 
Here 7r 7 (e) = V^P-y(e) and A 7 (e) = T ,(5) A 7 (e) are given 

by 

n a (e) = J de'p(e')f(e')Uz-e% (AT) 
n p (e) = J de'p(e')f(-e')^(e + e'), (A-5) 

Aa(e) = -C/ |A a (e+)| 2 7r a (e), (A-6) 
X fi (e) = -C/ 1 |A /3 (e+)| 2 7r /3 (e), (A-7) 
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which follow from eqs. (34a), (34b), (35a) and (35b), re- 
spectively. Equations (37) and (39) give expressions for 
We) =V<®R(, a (e) and f Q/3 (e) = V^R a p{e) as 

£p a (e)=Z f} (e)Re[A 2 a (e + )} 

+ 2Rc[i? /3 (e+)]Re[A ct (e+)]A a (e), (A-8) 

tap® = Ue)M4 

+ 2Rc[i? ct (e+)]Re[A /3 (e+)]A^(e). (A-9) 

Norms of the spectral functions of the defect propagator, 
namely the partition function Zf, are determined by the 
following sum-rule: 

fdeJ2^(e) = l- (A-10) 

J 7 

This equation follows from eq. (7). 

Appendix B: Equations for the Magnetic Sus- 
ceptibility 

The magnetic excitation and the static magnetic sus- 
ceptibility is given by eqs. (51) and (52) in the forms 
including the Boltzmann factor, respectively. For numer- 
ical calculations, we rewrite eq. (54) as well as eqs. (51) 
and (52), following ref. 29. We introduce, for each value 
of uj, simplified notations as follows: 



p(e 
p(e 
r(e 
r(e 



= n M (e + ,e + ^ + ), 
= n M (e~,£ + w + ), 
= R 1 (e+)R 1 (e + w+), 
= R 1 (e-)R 1 (e + uj + ), 

= i? i(e + )Ai(e + uj+) + Roi(e + w+)A 1 (e+), 

= fl i(e")Ai(e + w+) + R i(e + w+)Ai( e -). (B-l) 



The argument to has been omitted in the left-hand side 
of each definitions. Equation (54) is rewritten with use 
of the above symbols as 

p{e)=r{e){l-V 2 U 1 J de" p(e")f(-e")q(e e") 

x / de'p(e')f(e')p(e + e'-e")\. (B-2) 



We further introduce an operator V defined by 



\P(e)-P(e)]- 



(B-3) 



With use of the operator V, the imaginary part of the 
magnetic susceptibility is given by 

Im XM (u + ) = JViC(l - e"^) J deReVp(e). (B-4) 



We derive an integral equation for the Vp(e). Operating 
V on eq. (B-2), we obtain 

Vp(e) = p(e)[Vr(e)]/r(e) 

- V 2 Uif(e) J dep(e")f(e")[Vq(e-e")} 

x J de'p(e')f(e')p(e + e'-e") 

- V 2 U\f{e) j dep(e")f(e")q(e-e") 

x J 'de'p(e')f(-e')Vp(e + e' - e"). (B-5) 

As an initial condition of this equation, Vr(e) and Vq(e) 
are evaluated from the original quantities as follows: 

Pr(e)=i?i(e)i?i(e + w + ), 

Vq(e) -i|oi(e)A 1 (e + o;+) + i J Roi(e + ^ + )Ai(e). (B-6) 

In the first step, p(e) is computed by the numerical itera- 
tion of eq. (B-2) for a certain value of w. Then the results 
are substituted into eq. (B-5) to compute Vp(e), which 
leads to the dynamical magnetic susceptibility through 
eq. (B-4). This procedure is iterated for each value of lo. 

For the static magnetic susceptibility, we introduce, 
instead of V, another operator V defined by 

Vp(e) = -^e^ e Imp(e). (B-7) 



TrZf 

Then xm (0) is evaluated with use of V as 
Xm(0) =NxC j deVp(e). 



(B-8) 



The equations for Vp(e) are obtained by replacing V, f 
and q with P, r* and q* , respectively, in eq. (B-5). An 
initial condition for the equation is given by 

Pr(e) = -2&(e)Refli(e), 

Vq(e) = -2| i(e)Ai(e+) - 2 J R* 1 (e+)A 1 (e). (B-9) 

Vp(e) is evaluated by combination of eq. (B-2) with w = 
and modified one of eq. (B-5). 
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